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Forest biomass is one of the renewable and sustainable sources of energy that can be used for producing 
electricity, heat, and biofuels. The complex supply chain of forest biomass for energy generation, which 
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consists of different players and products and is affected by biomass characteristics, such as low density 
and unpredictable quality, makes the energy generation cost from biomass higher than that of the 
conventional sources of energy, such as fossil fuels. Moreover, variability and uncertainty in this supply 
chain, mainly due to the nature of material, economic condition and market fluctuation, affect the 
Keywords: amount of produced energy and its cost. Mathematical modeling, in particular optimization techniques, 
Forest bioenergy supply chain _ can be employed to manage the supply chain and achieve the optimum design. This paper reviews 
Renewable and sustainable energies studies which used deterministic and stochastic mathematical models to optimize forest biomass supply 


e chains for electricity, heat and biofuels production. Optimization models were used to provide the 
optimum solution for decisions related to the network design, technology choice, plant size and location, 
storage location, mix of products and raw materials, logistics options, supply areas, and material flows. 
Mainly, economic objectives were considered in these models. Further studies should consider 
environmental and social objectives, in addition to the economic ones, in the models. In non- 
deterministic models uncertainty mainly in the demand, supply, prices, and conversion yields were 
incorporated. Although material quality is an important uncertain parameter in the forest biomass 
supply chain that affects the amount and cost of produced energy, its variation was not considered in 
previous studies. 

© 2013 Elsevier Ltd. All rights reserved. 
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1. Introduction 
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Increasing the contribution of biomass in energy generation is 
considered as an important step in developing sustainable com- 
munities and managing greenhouse gas emissions effectively [1]. 
Although the conversion and transportation of forest biomass for 
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Fig. 1. Forest biomass supply chain for energy production. 


energy generation affect the air quality negatively [2], energy 
generation from forest biomass has the potential to decrease 
carbon emissions significantly when it substitutes fossil fuels [3- 
6]. Converting forest biomass to energy also has the potential to 
recover the waste that would otherwise be disposed to landfills or 
be incinerated, create jobs, provide local and sustainable energy 
for communities, and decrease their dependency on the interna- 
tional fuel market. Compared to other renewable energy sources 
such as wind or solar, the advantage of using forest biomass for 
energy generation is that it can be stored and used on demand [2,7]. 

Despite these advantages, there are several barriers in utiliz- 
ing forest biomass as an energy source. One challenge is the 
complexity of forest bioenergy supply chains which is partly 
resulted from complexity in forest industries supply chain in 
general. Forest industries consist of different interrelated and 
interconnected sectors, products and markets [8]. As an example, 
the raw material of the bioenergy facilities can be supplied from 
final products or by-products of some other plants such as 
sawmills, and pulp and paper mills. Therefore, the economic 
situation and the production mix of these plants can directly 
impact the raw material availability and its quality for bioenergy 
production. Moreover, unlike fossil fuels, forest biomass is usually 
speared over large areas rather than being concentrated. Also, 
forest biomass is a bulky material with relatively low density 
(400 and 900 kg/m? [7]), and high moisture content [3]. These 
characteristics contribute to high cost and complexity of forest 
biomass supply logistics [7]. The transportation cost could 
account for 50% of the total delivery cost of biomass in some 
cases [9] and the logistics system could comprise large number of 
equipment pieces and different transportation modes [2]. It is 
usually needed to convert forest fuel into chips before delivering 
it to the customer [10]. Inaccessibility of forests in some months 
during the year, when energy demand is quite high, raises 
concerns about the secure supply of biomass to energy plants. 
Therefore, storage, which affects the quality of material [11], is 
also important in this supply chain. Comminuting and storing 
residues can be done either in the forest, at the plant or at an 
intermediate point. The location of an energy plant also plays a 
key role in the economic performance of energy generation from 
forest biomass. There is variability and uncertainty in forest 
bioenergy supply chains due to several other factors, such as 
market instability, natural disasters, and policy and climate 
changes as well as the nature of the industry (for example 
heterogeneous raw material and unpredictable quality [3]). 
Uncertainty makes this supply chain volatile and risk vulnerable, 
which in turn makes the proper planning difficult. All these 
challenges contribute to high cost of forest bioenergy compared 
to other sources of energy. Utilizing more advanced technologies, 
for example to improve raw material quality or system efficiency, 
is one way to deal with some of these challenges. Another 
solution for improving the performance of a forest bioenergy 
supply chain is to optimize its design and management. 


Mathematical programming models can be used to optimize 
these supply chains. These models are effective particularly when 
different parts of the supply chain, such as procurement, produc- 
tion, transportation, and distribution, and different decision levels, 
such as strategic, tactical and operational levels, are integrated. 
Usually, decisions regarding the supply chain design such as 
location, technology and capacity are made in the strategic level, 
while decisions related to the flow of material and production 
planning can be made in tactical and operational levels. 

Optimization techniques have been employed for modeling supply 
chains in different industries including forest industries [8]. In addition 
to modeling the supply chain, optimization has the advantage of 
providing the optimum solution based on the objective function 
(s) defined in the model. Using optimization techniques in designing 
and management of forest bioenergy supply chain can result in better 
performance which helps making this energy source economically 
viable [12]. In several previous studies, optimization techniques have 
been employed to manage the forest bioenergy supply chain for heat, 
electricity and biofuels production from strategic, tactical and opera- 
tional point of views. Most of these studies were deterministic and 
ignored uncertainty, while there are examples that included uncer- 
tainty in the supply chain models especially during the past few years. 
This paper reviews all of these studies. 

A number of papers reviewed the literature related to different 
parts of the bioenergy supply chain. Jebaraj and Iniyan [13] 
reviewed energy models and allocated a small section of renew- 
able energies in their paper to biomass and a section to optimiza- 
tion models in energy systems in general. Bafios et al. [14] 
reviewed studies that used optimization methods in renewable 
and sustainable energies and dedicated a section to the bioenergy 
industry. These studies were not comprehensive and only pro- 
vided examples from the bioenergy supply chain. Wang et al. [15] 
and Scott et al. [16] reviewed the multi-criteria decision making 
methods applied in sustainable energy decision-making and 
bioenergy systems, respectively. These studies did not focus on 
optimization and supply chain design and management problems. 
Johnson et al. [17] studied methods and literature on optimum 
location of forest biomass—biofuel facilities and did not consider 
other decisions in the supply chain. The potential application of 
geographic information system (GIS) in evaluating the feasibility of 
bioenergy projects was studied by Calvert [18]. As noted by the 
author, geographical aspects are important factors affecting the 
feasibility of such projects. The use of spatial data and GIS in 
modeling of the value chains could be helpful in location analysis, 
transportation cost estimation, and actual (vs. potential) feedstock 
availability quantification. Awudu and Zhang [19] reviewed studies 
that considered uncertainty in biofuel supply chains. The authors 
only focused on the biofuel industry and did not include determi- 
nistic cases or those related to heat and energy plants. Moreover, 
they considered both agricultural biomass and forest biomass 
which have different origin and different supply chains. The 
present paper is more comprehensive compared to the above 
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mentioned review papers since it covers previous studies on 
optimization of forest bioenergy supply chain design and manage- 
ment, including those for heat, electricity and biofuels production. 
It also discusses the issue of uncertainty in these supply chains, its 
sources and the methods used for dealing with it. This paper gives 
a comprehensive review on the deterministic and stochastic 
models used for forest bioenergy supply chain in literature. The 
studies are categorized in two groups: (1) studies that used 
deterministic mathematical programming for modeling biomass 
supply chain in generating heat and/or electricity (mostly in 
district energy systems) and biofuels, and (2) studies that incor- 
porated uncertainty in their modeling of forest bioenergy supply 
chain. The main focus of the reviewed studies was on optimizing 
the supply chain profit, while other objectives related to environ- 
mental and social issues should be considered in future studies. 


2. Deterministic optimization models 


The wood flow in the forest products supply chain starts from 
harvest areas to value added mills such as sawmills, pulp and 
paper mills, secondary wood products, wood pellet mills, and 
eventually bioenergy plants [8]. These mills produce different 
products, and may generate co-products (demand-driven) or by- 
products (secondary results of production). Forest biomass for the 
purpose of energy generation can be supplied from forest residues 
including branches and tops left in the harvest areas, by-products 
of other forest product mills, such as sawdust, bark and shavings 
[7], and fast growing crops such as poplar and willow grown 
specifically for energy purposes [20]. Chipping, handling, trans- 
porting, storing and pre-processing operations, such as drying for 
improving the quality of biomass, are usually needed before using 
the forest biomass for energy generation. The energy production 
process depends on the conversion technology used in the energy 
plant [21]. Forest biomass can be used for energy generation either 
directly, such as in direct heat and power generation, or indirectly 
such as in biofuels (pellet, bioethanol) production. The energy 
generation technologies include pelletization, combustion, co- 
combustion, gasification, pyrolysis, digestion and fermentation 
[7,22]. Different types of energy products are then sent to 
customers through the grid, networks or channels of distributors, 
wholesalers and retailers. Profit in each stage of the forest 
bioenergy supply chain is a function of procurement, transporta- 
tion, operating, capital and other costs, and depends also on the 
availability and quality of biomass [23]. Fig. 1 shows the network of 
generating energy from forest biomass. 

Different optimization techniques, such as linear programming 
(LP) and mixed integer linear programming (MILP), have been 
used for supply chain design and management. LP is a mathema- 
tical method which includes a set of variables to be determined, a 
linear objective function to be optimized, and a set of linear 
equality or inequality constraints to be met. The main advantages 
of using this optimization method are the ability of LP algorithm to 
solve large scale problems, its assured convergence to global 
optimum solutions, having no need to have an initial solution 
and its use of a well-developed duality theory for sensitivity 
analysis and the ease of problem formulation. If some of the 
variables in LP are integer, the model is called mixed integer linear 
programming (MILP). In this section, the studies that used opti- 
mization in supply chain design and management of district 
heating systems, electricity plants, biofuel plants and co- 
generation plants that utilize forest biomass are reviewed. 


2.1. District heating plants 


One of the main uses of forest biomass for energy generation is 
in district heating systems. In a district heating system (DHS), a 


central plant generates thermal energy and a network of pipelines 
distributes the produced energy in the form of hot water or steam 
to a group of customers in a community [24,25]. Heat exchangers 
installed at each building draw off the heat from water and when 
heat has been circulated within the interior atmosphere, a sepa- 
rate pipe network returns the water to the central plant for 
reheating [26]. In addition to heat and hot water, district energy 
systems may provide chilled water. 

District heating systems became popular after the oil crisis in 
the 1970s because of their high efficiency in energy generation 
[26] compared to decentralized systems, which produce and 
consume energy within the same building. As a consequence of 
higher efficiency, primary fuel consumption decreases and some 
negative environmental impacts of burning fuels might be offset. 
District heating systems can use a wide variety of fuels including 
fossil fuels and renewables and offer better pollution control than 
decentralized systems [27]. 

Some studies have shown that the costs of using biomass for 
energy generation in district energy plants are higher than other 
sources of energy such as fossil fuels [28-32]. While this fact may 
seem to be discouraging for making investments in biomass 
district heating projects, it is notable to consider that forest 
biomass utilization for power and thermal energy generation has 
environmental, social, and economic benefits in some cases. 
Governmental subsides and grants could help these plants to be 
competitive with other sources as was the case in Sweden, a 
pioneer country in exploiting forest biomass in energy plants [33]. 
The effect of biomass price and/or subsidies on the profitability 
and structure of heating plant supply chains have been studied 
and discussed in [34-37]. 

The cost of energy generation in a district heating plant is 
related to the delivery cost of biomass, location of plants, the type 
of technology that is used to convert biomass into a useful from of 
energy and the operating scale of energy generation plants. The 
cost advantage of larger systems, known as economy of scale, is a 
known fact. However, larger systems require more forest fuel and 
consequently more raw material which result in having higher 
transportation costs. A significant cost factor in energy generation 
from forest biomass is the conversion technology cost. Purchase 
price, installation, pre-treatment, operating and maintenance, and 
labor costs differ for various types of technical solutions. Further- 
more, conversion technologies are different in terms of efficiency 
which dictates the amount of required biomass and affects the 
transportation cost. 

Choosing the conversion technology is one of the strategic 
decisions that may be made in the design of a heating plant supply 
chain. For example, decisions can be made to select between a co- 
generation system and an individual system [34], or select among 
different technologies, such as grate firing combustion (GFC), 
fluidized bed combustion (FBC), fluidized bed gasification (FBG), 
and fast pyrolysis (FP) [22]. In [34], the investment cost, fixed and 
variable costs, fuel cost, waste disposal cost and revenue from 
selling energy (heat and/or electricity) were the components of 
the objective function of a mixed integer programming model. In 
[22], a linear programming optimization model was developed for 
each individual technology and the optimum values of their 
objective functions were compared. The authors concluded that 
despite the low efficiency of grate firing combustion, it was the 
best technology to be adopted because of its lower installation and 
Management costs compared to other alternative technologies. 
The same authors then changed their model to mixed integer 
programming to find the optimal conversion technology with only 
one run of model using binary variables in the model [38]. 

Usually, an optimization model can be used to make decisions 
about different variables, such as size and location of the plant, 
material flow, etc. at the same time. Chinese et al. [35] used a 


Table 1 


Summary of studies on deterministic optimization of forest biomass district heating plant supply chains. 


Author 


Eriksson et al. 
[39] 


Nagel [34] 


Gunnarsson et al. 
[40] 


Chinese et al. [35] 


Frombo et al. [22] 


Frombo et al. [38] 


Kanzian et al. [42] 


Van Dykenetal. 
[43] 


Keirstead et al. 
[41] 


Han and 
Murphy [44] 


Year—region Objective function 


1989—Sweden 
costs) 


2000—Germany Maximizing annual profit (revenue from sale of energy minus investment cost, 


fixed and 
variable costs, fuel cost and waste disposal cost) 


2004—Sweden 


2005-—Italy Maximizing annual profit (revenues from sale of thermal energy and charging 
customers with connection fees minus boiler investment cost and construction 


and operating costs of facility) 


2009-—Italy Maximizing net annual profit (revenue from sale of heat and power minus 
felling and processing, skidding, highway transportation, plant installation and 


management costs) 


2009-—Italy Maximizing net annual profit (revenue from sale of heat and power minus 
felling and processing, skidding, highway transportation, plant installation and 


management costs) 


2009—Austria 
transporting costs) 


2010—Norway 
salvage value) 


2012—UK Minimizing system cost (biomass purchase, storage, transportation and 


conversion costs) 


2012—US Minimize the weighted sum of transportation costs 


Minimize the total working time 


Minimizing supply cost of forest biomass (chipping, storing and transportation 


Minimizing biomass supply cost (transportation, chipping and storage costs) 


Minimizing biomass supply cost to the heating plants (chipping, storing and 


Minimizing the present value of the costs (investment and operating costs and 


Decision variables 


Flow of biomass direct or via storage Chipping location 


Level of heat produced by each boiler at each time period 

The capacity of the system 

Whether or not to integrate a boiler into the heating system (Binary 
variable) 


Flow of biomass within the supply network 

Quantity of biomass chipped and stored at roadside and terminal 
If biomass is forwarded to or is chipped at each roadside location 
(Binary variable) 

If each sawmill is contracted or not (Binary variable) 

If each terminal is used or not (Binary variable) 


Heat produced by each boiler at each time period 

The capacity of the system 

If a boiler would integrate to the heating system or not (Binary 
variable) 


Annual quantity of biomass harvested from each supply area 
The plant capacity for different conversion technologies 


The quantity of biomass harvested at each harvesting location and to be 


used at each plant location 
The capacity of each plant 
Selection of the conversion technology (Binary variables) 


Volume of wood chips transported from each terminal to each plant 
Location of terminals and plants (Binary variable) 


Biomass and product flow within the supply network 

Energy flow within the supply network 

Emissions from storing and drying biomass 

Biomass input and output moisture content to and from dryer 
Biomass storage duration (Binary variable) 


Optimal capacity of boilers 

Whether chipped forest biomass should be imported from neighbor 
area or non-chipped residues should be imported and then chipped 
within the area (Binary variable) 


Truck schedules 


Method 


Linear Programming (LP) 


Mixed Integer 
Programming (MIP) 


Mixed Integer 
Programming (MIP) 


Mixed Integer 
Programming (MIP) 


Linear Programming (LP) 


Mixed Integer 
Programming (MIP) 


Mixed Linear 
Programming (MIP) 


Linear and Mixed Integer 
Programming (LP and 
MIP) 


Mixed Integer 
Programming (MIP) 


Simulated annealing 


COE 
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mixed integer model to obtain the most profitable configuration of 
a multi-source biomass district heating plant including decisions 
about the optimal size of the plant, resources to be integrated to 
the network and the sequence of integrating boilers to the net- 
work under different government subsidy levels. The energy 
sources consisted of the exhausted steam from an existing district 
electricity generation turbine, along with wood wastes from local 
manufacturers. Five individual biomass-fired boilers with different 
capacities were identified to provide process heat for this indus- 
trial district. Virtual boilers including a new natural gas boiler and 
a new woody biomass boiler were also considered. 

Finding the optimal flow of biomass from buffer storages to 
energy plants while minimizing the supply cost has been of great 
interest in the available literature. At the tactical level, one of the 
premier studies which used mathematical programming in the 
operation of a district heating system was done by Eriksson and 
Bjérheden [39]. They developed a model to make the following 
decisions: whether to transport biomass to the heating plant 
directly or via storage, and whether to chip biomass at roadside 
or at the plant. The results of the optimal solutions were to chip 
residues at the roadside and direct flow of chips to the heating 
plant. When storage was used, additional transportation cost 
incurred and the authors showed that this additional transportation 
cost would not be paid off by the improved quality (better moisture 
content) of biomass after storage. A more complex situation where 
several storages and heating plants are included in the supply chain 
could be dealt with decision support systems (DSS) as was imple- 
mented in [40]. The objective of this mixed integer based DSS was to 
minimize the supply cost of biomass to heating plants including 
forwarding, chipping, storing and transportation costs. The DSS 
defined the optimal location for chipping forest residues, and an 
optimal monthly plan for flow of wood chips from roadside to the 
storages and plants, from storages to heating plants, and also from 
sawmills to the storages and plants. The mixed integer programming 
model in a study [41] found the optimum cost of forest biomass 
supply for fulfilling the heat demand of an urban area. This 
optimization model was used to decide between importing wood 
chips and forest residues from nearby area. When the decision was to 
import forest residue, the chipping process was carried out within 
the studied area. The objective function in this study was to minimize 
the annual cost including annual fuel, biomass storage, transporta- 
tion, and conversion costs. A supply chain of 16 combined heat and 
power plants and 8 terminal storages with a total capacity of 
178,250 m? in Austria were studied in [42]. The authors presented 
an optimization model to minimize the biomass supply cost to the 
plants, which included chipping, storing, and transporting costs, 
while defining the optimal flow of wood chips. The results of this 
study were similar to those in [39], which indicated direct flow of 
biomass was less expensive than indirect flow via terminal. 

Van Dyken et al. [43] developed a linear mixed-integer model for 
biomass supply chain with transport, storage and processing opera- 
tions at the operational level. In their case study, they considered 12 
weekly time steps and three different biomass products including 
spruce, wood chips and wood pellets. There were constant weekly 
demands for both wood chips and heat. The demand for wood chips 
was fulfilled partly from chip supply and partly from chipping spruce 
in a chipper with an additional cost. To improve the biomass quality, 
two kinds of drying were considered: active drying through dryer and 
passive drying through storage. Active drying increased the cost and 
emissions since it consumed energy, while passive drying caused 
volume loss. After drying, part of chips was used to fulfill biomass 
demand and the rest was sent to the second storage which followed 
into a pellet mill and a heating plant. The optimization model was 
solved and the optimum operational plan, which minimized the total 
cost including the fuel cost, biomass cost and electricity cost, was 
derived for the case study. 


A truck scheduling optimization problem was studied in [44] 
for transporting forest biomass from harvesting sites and sawmills 
to heating plants, pulp plants or export harbors in Oregon, US. The 
problem contained four types of forest biomass including hog fuel, 
chips, sawdust, and shavings transported from 45 sawmills to 20 
conversion facilities over a 1 year period using 75 chip vans and an 
average of 100 truckloads per day. The objective function was to 
minimize the weighted sum of transportation costs and the total 
working time. Different constraints were included in the model 
such as routes, working time and predetermined order require- 
ments. A heuristic algorithm, called simulated annealing, was used 
for solving this optimization problem. The author indicated that 
this algorithm provides “near optimal” solution which was still 
better than the existing situation in terms of transportation cost 
and total travel time (18% and 15% lower, respectively). 

Table 1 summarizes the studies discussed in this section. 


2.2. Power plants 


Forest biomass can be used in power plants for generating 
electricity. It can be burnt at a constant rate in a boiler furnace to 
heat water and produce steam. Then, the steam is carried through 
the furnace using pipes to raise its temperature and pressure 
further. Finally, the steam passes through the multiple blades of a 
turbine, spinning the shaft and the shaft runs an electricity 
generator which produces an alternating current to use locally or 
to supply the national grid [7]. 

The optimal supply area and location of a forest biomass power 
plant in a distributed power generation system was determined by 
Reche et al. [45]. The objective function was to maximize the 
profitability index as a function of the net present value of benefits 
from the sale of electrical energy minus the initial investment, 
collection, transportation, maintenance and operating costs. The 
authors used an artificial intelligence method, called particle 
swarm optimization. They concluded it is important in distributed 
generation systems to consider the technical constraints of the 
network and the voltage regulation. Finally, they evaluated the 
model performance using simulation. 

Alam et al. [46] constructed a three-objective model for optimizing 
the amount of each individual type of biomass from each of the 
harvesting zones, and then applied their model to a 50MWh 
biomass power plant using both harvesting residues and poplar trees 
collected from three management zones in Northwestern Ontario, 
Canada. To optimize the supply chain of energy plants, it is sometimes 
necessary to formulate a problem with more than one objective since 
single objective models cannot always represent the problem accu- 
rately. The objectives are often in conflict (minimizing and maximiz- 
ing objectives) and it might not be possible to achieve an optimal 
solution that optimizes all the objectives simultaneously. In this 
situation, the trade-off between objectives can be shown and the 
most efficient solution is selected. In [46], pre-emptive goal program- 
ming was applied to give priorities to the objectives as follows: (1) 
minimizing the procurement cost of feedstock, (2) minimizing the 
transportation distance of biomass supply to the plant, and (3) 
minimizing the feedstock moisture content. 

Finding the optimal size, location, supply area and net present 
value of an electricity plant in Spain was studied in [47]. The raw 
material of the power plant was olive tree pruning residues and 
the technology for electricity generation was gasifier with gas 
turbine. The authors used GIS data for the location and number of 
olive trees per kilometer square, roads, topographical features, 
electric lines location, etc. Different plant sizes and locations were 
considered and the optimal one with the highest net present value 
was determined using three metaheuristic methods. These meth- 
ods were Genetic Algorithms (GA), Binary Honey Bee Foraging 
(BHBF) and Binary Particle Swarm Optimization (BPSO). It was 
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Table 2 
Summary of studies on deterministic optimization of forest biomass power plant supply chains. 
Author Year—region Objective function Decision variables Method 
Reche et al. 2008—Spain Maximizing profitability index (net present value of revenue from selling Location and supply area of the Particle swarm 
[45] electricity minus initial investment, biomass collection and transportation costs, biomass power plant optimization 
and maintenance and operation costs) 
Alam et al. 2010—Canada Minimizing total biomass procurement cost Quantity of biomass procured from Multi- 
[46] Minimizing total distance for procurement of biomass each supply location to each plant objective 


Maximizing the quality of biomass (minimizing moisture content) 


Vera et al. [47] 2010—Spain 


costs) 


Shabani and 2013—Canada 


Maximizing net present value (revenue from the sale of electrical energy minus Pant size and location 
initial investment and collection, transportation, maintenance and operation 


Maximize the profit (revenue from selling electricity minus procurement, 


Biomass procurement zone 
selection (Binary variable) 


Programming 


Several 
metaheuristic 
methods 


Supply area 


How much biomass to purchase, Mixed Integer 


Sowlati [48] transportation, production, storage and ash management costs) store and consume in each month Programming 
Whether or not to produce (MIP) 
surplus electricity (binary) 

Table 3 
Summary of studies on deterministic optimization of forest biomass biofuel plant supply chains. 
Author Year—region Objective function Decision variables Method 
Chinese 2009-—Italy Minimizing total cost of supply chain Flow of biomass within the supply network Mixed Integer 
et al. [51] (harvesting, transportation, processing Whether to use a pre-processing equipment, Programming 
and facility installation costs) e.g. dryer, or not (Binary variable) (MIP) 
Eksioglu 2009—USA Minimizing total annual cost (investment, Number, size, and location of bio- refineries required Mixed Integer 
et al. [52] harvesting, storing and transportation costs) Quantity of biomass harvested, shipped, processed and stored Programming 


Whether a biorefinery and a collection facility with specific size (MIP) 
are located in each site (Binary variables) 


Santibañez 2011—Mexico Maximizing profit (revenue from sale of 
et al. [53] products minus investment, process, operating 
and transportation costs) Minimizing 
environmental impacts 


concluded that BHBF algorithm converged to the optimal solution 
better than BPSO and GA. The results indicated that the optimal 
plant size was 2 MW and the predicted optimal location of the 
plant was in the area with highest available biomass. 

In a recent study, the value chain of a forest biomass power 
plant was optimized [48]. The authors included forest biomass 
procurement, storage, production and ash management in a 
dynamic model consisted of a 1-year plan with monthly time 
steps. The objective function was to maximize the power plant's 
profit, while meeting the constraints related to production capa- 
city, maximum available supply, demand fulfillment, material flow 
balance and storage limits considering biomass quality, in terms of 
moisture content and higher heating value, in the optimization 
model. The problem was solved using a mixed integer non-linear 
model and was applied to a real case study. Different scenarios on 
biomass availability and investment in a new ash recovery system 
were considered. It was concluded that investment in a new ash 
recovery system had economic and environmental benefits for the 
power plant. The sensitivity analysis results showed that variations 
in some parameters, such as maximum available biomass, elec- 
tricity price, biomass price and quality, had higher impact than 
other parameters on the optimum solution. 

Table 2 summarizes the studies using optimization models in 
forest biomass power plant supply chains. 


2.3. Biofuel plants 
The bioethanol production has increased in recent years in many 


countries, such as US. Although most of the bioethanol is produced 
from agricultural biomass, the controversial issue of using plants as 


The quantity of products produced from different biomass 
feedstock using different processing routes 

The quantity of each biomass feedstock used for producing 
different products through different processing routes 


Multi-objective 
Programming 


fuel instead of food made it necessary to look for more acceptable 
sources, namely forest biomass [49]. Generating biofuels from forest 
biomass is still in the developing phase and has not been commer- 
cialized yet. The main challenges in commercialization of this 
technology include high energy consumption for woody biomass 
pretreatment, even when compared to agricultural biomass, low 
system efficiency, process scalability and intensive capital investment 
[50]. In most of the studies presented here, forest biomass combined 
with agricultural biomass is used for biofuel production. 

Chinese et al. [51] considered a real-life problem of supplying a 
biofuel plant with forest fuel. A mixed-integer linear programming 
model was proposed to determine the optimal configuration of the 
supply chain. It was mentioned that the model could be useful in 
resolving trade-offs between decentralized early treatment of 
biofuels, resulting in lower transportation costs, and centralized 
final treatment, allowing to reap the benefits of economies of scale. 
It was therefore advised to apply integrated supply chain planning 
concepts to design biofuel logistics systems and to support policy 
making in the energy field. An MIP model was also developed by 
Eksioglu et al. [52] for designing the bio-refinery supply chain 
producing cellulosic ethanol from agricultural and woody biomass. 
The model outputs were the number, size and location of bio- 
refinery plants with the objective of minimizing the total cost of 
annual harvesting, storing, transporting and processing biomass; 
storing and transporting ethanol; and locating and operating bio- 
refineries. The model included constraints on biomass availability, 
flow, conversion, production and inventory capacities, and 
demand. The data from the State of Mississippi was used to 
validate the model. The author concluded that transportation 
costs, biomass availability, technology type, and planting and 
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harvesting costs are important factors in supply chain design 
decisions. 

The trade-off between economic and environmental objectives 
in the optimal planning of a bio-refinery in Mexico was evaluated 
in [53]. The authors used a multi-objective optimization model for 
selecting the feedstock type, processing technology, and a set of 
products in a bio-refinery supply chain. The raw material con- 
tained different types of agricultural biomass, wood chips, saw- 
dust, commercial wood for producing ethanol, hydrogen, and 
biodiesel (generated only from agricultural biomass). The objec- 
tives were (1) to maximize the profit considering the costs of 
feedstock, products, and processing, and (2) to minimize the life 
cycle environmental impacts. The authors applied their model to a 
case study in Mexico. The decision makers could select from the 
output the solutions that fit the specific requirements and com- 
pensate for both objectives simultaneously. 

Table 3 summarizes the optimization studies in forest biofuel 
plant supply chains. 


2.4. Co-generation plants 


Combined heat and power (CHP) systems are another type of 
energy plants that use co-generation technology to produce both 
heat and power in a district heating system [25]. CHP systems have 
up to 35-40% higher efficiency than conventional power genera- 
tion systems [27]. In some studies, mathematical programming 
was used to compare the cost of generating either energy 
or biofuels from biomass and also to evaluate the possibility of 
co-generation. 

Some studies indicated that utilizing biomass for energy gen- 
eration is more cost efficient than for biofuel production. Azar et al. 
[54] proved that utilizing biomass in the heat sector was the most 
economical scenario. Wahlund et al. [55] showed that using wood 
biomass for pelletization would have a lower cost and higher CO3 
reduction than using it for biofuel production. Feng et al. [56] 
investigated the possibility of having bioenergy facilities (they 
called it biorefinery) in typical existing sawmills, pulp and paper 
mills, wood panel facilities, biochemical, energy, and pellet facil- 
ities. Then, the authors developed a mathematical model to design 
this integrated supply chain optimally. 

A methodology for optimizing the utilization of distributed 
biomass resources for energy production was proposed by Alfonso 
et al. [57]. The main focus of the paper was to optimize the 
logistics, but economic and environmental analyses of different 
bioenergy alternatives were also performed. The authors indicated 
that the methodology would provide the optimal locations of the 
biomass plant, energy application (electricity, heat and/or stan- 
dardized biofuels such as pellets), and the employed technology. 
This methodology was applied to three districts in Spain. Based on 
the results, the authors concluded that the shortest payback period 
and highest CO, savings were attained from co-generation plants, 
followed by pellet plants. The least ranked option was power-only 
power plants. 

Difs et al. [36] analyzed different biomass gasification scenarios, 
and determined the optimum configuration with the current fossil 
fuel price and green energy policies. Wetterlund and Söderström 
[58] considered two scenarios of co-generating Synthetic Natural 
Gas (SNG) and district heat, and co-generation of heat and power. 
The authors determined the policy support levels (tradable biofuel 
certificates) that would make the SNG scenario cost competitive 
with CHP, while maximizing the annual profit over a 20-year time 
period. They concluded that no financial support would be 
required for combined heat and power generation because the 
market price of electricity was high enough. 

The cost-effectiveness of different applications of biomass 
gasification was analyzed by Börjesson and Ahlgren [37]. The focus 


of this study was to determine whether CHP generation in biomass 
integrated gasification combined cycle (BIGCC) plants, and biofuels 
production in biomass gasification biorefineries in a case study in 
Sweden were cost efficient. Examples of the model outputs 
included electricity production, use of electricity, biofuels produc- 
tion, required biofuel subsidy level, CO2 emissions, and total use of 
biomass. The authors indicated that the outcome of their model 
was not an optimized system because the systems they studied 
consisted of a number of subsystems, each controlled by different 
actors with probably a bigger interest of optimizing their own 
subsystem than the overall system. They mentioned that the 
results were rather an illustration of what could be than a 
prediction of the future development. 

Schmidt et al. [59] used a mixed integer linear programming 
model for optimizing the location of bioenergy plants using forest 
biomass in Austria. The bioenergy plants included integrated 
gasification combined cycle (IGCC) system and biomass CHP plants 
with carbon capture storage (CCS), pellet mill, and transportation 
fuel (methanol and ethanol) plants. The optimal size of the plant 
and biomass allocation area for different systems were also 
considered. It was found that the production of transportation 
fuels were not economical. 

The problem of indicating whether to produce electricity in 
addition to heat at biomass combustion plants was studied by 
Freppaz et al. [60]. Their model's objective function was to 
maximize the financial yield of six proposed plants, while deciding 
about the annual amount of biomass harvested in each collection 
area, the annual flow of biomass from each collection area to each 
plant, the percentages of heat and electricity generated in each 
plant and the capacity of combustion system in each plant. 

A decision support system (DSS) was presented by Rentizelas 
et al. [61] to optimize a multi-biomass energy conversion system 
to generate electricity, heating and cooling in an area in Greece. 
The objective function was to maximize the financial yield of the 
investment, while meeting customer demand. The decision vari- 
ables included the location and size of the bioenergy facility, the 
biomass types and quantities and the maximum collection dis- 
tance for each type. The authors concluded that considering multi- 
biomass supply chain reduced the cost by decreasing warehousing 
requirements, especially for seasonal types of biomass. The devel- 
oped model was non-linear and a hybrid optimization method was 
used to solve that. 

Rauch and Gronalt [62] developed a model for designing a 
forest fuel CHP plant supply chain in Austria by making decisions 
about transportation modes and spatial arrangement of terminals. 
The model was a mixed integer programming model with the 
objective function of minimizing the total procurement cost. The 
effect of changes in forest fuel supply (domestic resources or 
imports), transport modes (truck only, truck and ship, or truck, 
ship and rail), energy price (increase by 0%, 20% and 40%), and 
truck load capacity (50%, 40% and 30%) on the overall cost was 
evaluated. Eight scenarios were constructed and compared. Table 1 
shows the summary of studies discussed above. 

The summary of studies on optimization of forest biomass 
co-generation plant supply chains is provided in Table 4. 


3. Uncertainty in the forest biomass supply chains 


Uncertainty is referred to lack of information or lack of degree 
of belief in the validity of the information about the existing or 
future state of a system [63]. Uncertainty can result from measure- 
ment errors and ignorance, which is to some extent inevitable and 
might be reduced by further studies or investing in improved 
technology to acquire high quality data [64,65]. Uncertainty may 
result from variability in random future events due to their 


Table 4 


Summary of studies on deterministic optimization of forest biomass co-generation plant supply chains. 


Author 


Freppaz et al. [60] 


Alfonso et al. [57] 


Rentizelas et al. [61] 


Börjesson et al. [37] 


Difs et al. [36] 


Schmidt et al. [59] 


Wetterlund and 
Söderström [58] 


Rauch and Gronalt [62] 


Year—region 


2004—Italy 


2009—Spain 


2009—Greece 


2010—Sweden 


2010—Sweden 


2010—Austria 


2010—Sweden 


2011—Austria 


Objective function 


Maximizing annual profit (revenues from sale of energy 
minus harvesting, transportation, installation and maintenance, and energy 
distribution costs) 


Minimize transport duration, optimize 
the location, etc. 


Maximizing the financial yield of the investment 


Not discussed in the paper 


Maximizing annual profit (revenues from 
sale of energy products minus investment, 
fuel and maintenance costs) 


Minimizing total cost of energy generation 

(biomass supply and transportation costs, bioenergy generation cost, carbon 
capture and storage cost, 

plant building and distribution network investment costs, and distribution cost) 


Maximizing annual profit (revenues from sale of electricity and 
synthetic natural gas minus investment, fuel and maintenance costs) 


Minimizing total procurement cost (transport, chipping investment, operations 
and maintenance costs) 


Decision variables 


Annual quantity of biomass harvested at each 


collection area and transported from each collection area to each of 


six district energy systems 

Capacity of each plant 

The percentage of thermal energy generated at each plants 
If the plant produces electricity or not (Binary variable) 


Biomass resources, logistics structure, bioenergy plants size and 
location, technology type, etc. 


Location and size of the bioenergy facility 
The biomass types and quantities 
The maximum collection distance for each type 


The optimal production capacity at different subsidy levels 
Selection of alternative technologies for district heat generation 
(Binary variable) 


Capacity of new investment 
Selection of investment alternatives for future (Binary variable) 


The annual amount of energy commodities produced at plants: 
heat, power, pellets, and transportation fuels 

Plant size and location (Binary variable) 

Pipeline networks selection (Binary variable) 

District heating networks selection (Binary variable) 


The optimal government support level (subsidy) 
Selection of new investment alternatives (Binary variable) 


The annual volume of fuel transported between districts, 
terminals, regional departure railway, and the CHP plant 
Open or close a terminal (Binary variable) 


Method 


Mixed Integer 
Programming (MIP) 


Did not mention 


Hybrid optimization 


Mixed Integer 
Programming (MIP) 


Mixed Integer 
Programming (MIP) 


Mixed Integer 
Programming (MIP) 


Mixed Integer 
Programming (MIP) 


Mixed Integer 
Programming (MIP) 


90€ 
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inherent nature (such as feedstock characteristics) [66], which can 
be controlled to some extend by employing better forecasting 
methods and/or using expert judgment. It can also result from lack 
of reliable historical data or lack of certainty in historical data, for 
example lack of data on the demand of a new product. Other 
sources of uncertainty include imprecision in judgment, vague- 
ness, and ambiguity related to the known objects, which belong to 
poorly defined sets so they cannot be classified well [63-65]. This 
form of uncertainty is secondary as it may be the consequence of a 
lack of control in the manufacturing system, which can be 
decreased using some technologies such as part-tracking systems. 

From the system boundary point of view, the source of 
uncertainty may exist outside the production process, called 
environmental uncertainty, such as uncertainty in demand and 
supply, or it may be within the production process, called system 
uncertainty, such as uncertainty in lead time due to machine 
failure [67]. 

In terms of time horizon, uncertainty may be the result of short 
term variations, such as day-to-day processing variations, can- 
celed/rushed orders and equipment failure, or long-term varia- 
tions, such as raw material/final product unit price fluctuations, 
seasonal demand variations and technology changes. Therefore, 
uncertainty exists in supply chains at strategic, tactical and 
operational levels and should be considered in supply chain 
decisions. 

While some of the sources of uncertainty in forest biomass 
supply chains are similar to other industries, such as economic 
fluctuation and instability, raw material supplies, manufacturing 
process time, machine breakdown, reliability of transportation 
channels, and exchange rates, there are other sources of uncer- 
tainty that are related to the specific characteristics of forest 
biomass supply chains which are summarized here. 

Interdependency between different forest sectors: There are 
interdependencies between different sectors and markets within 
the forest industry supply chains. Consequently, variations in one 
part of the supply chain usually propagate into the other parts. For 
instance, the downfall of the economy and the US housing market 
collapse in 2007 resulted in the closure of many Canadian sawmills 
in 2008 [68]. This impact was not limited to the lumber industry 
and resulted in the decrease of residual chip supplies for pulp and 
paper mills, bioenergy power plants and other wood chip users. 
Some bioenergy power plants had to curtail their operation for 
several months due to the shortage of supply. 

Changing feedstock supply: The need for having a continuous 
supply of raw material for a bioenergy facility necessitates the use 
of a mixture of material or even to have new material sources in 
the future. Therefore, this industry has a dynamic supply chain 
over time which may have significant variations, such as variations 
in its environmental impacts [66]. 

Wood is a heterogeneous natural material: Its physical and 
chemical characteristics affect products quality and quantity [66]. 
In bioenergy industry, the moisture content and heating value of 
raw material play an important role in the amount of produced 
energy and its costs [6]. Heating value and moisture content vary 
from one tree stands or species to another [7,69,70] and also 
differs in different types of biomass (e.g. bark, sawdust, shavings) 
[71]. Wood properties may be affected by external factors, such as 
growth condition, climate, harvesting methods, storage and trans- 
portation methods. Biomass quality, such as moisture content, can 
also change during storage, production, and transportation. 

Divergent production structure: Forest products industries 
generally have a divergent production structure, which means 
multiple products, by-products, and co-products are made from a 
single product simultaneously. Consequently, it is difficult to 
completely control the manufacturing processes. Moreover, it is 
challenging to forecast the quality and quantity of outputs due to 


this production structure and the use of heterogeneous natural 
raw material in the production. This fact can impact the amount of 
raw material available for bioenergy plants, which are supplied 
from other forest product mills. 

Ambiguous values and objectives: Most forest areas include 
large areas with diverse geographical and ecological characteris- 
tics. In forestry, it is usually needed to incorporate different values 
and stakeholders’ preferences and interests which sometimes 
cannot be understood, interpreted and quantified completely 
[65]. Therefore, it is likely to have vague factors, values and 
objectives which can also exist in the forest bioenergy supply 
chain. This aspect of uncertainty cannot be dealt with like other 
sources of uncertainty. To some extent, it is possible to spend time 
and money in some forms of consulting with the stakeholders to 
get a better understanding of their preferences, opinions, and 
values. However, sometimes the stakeholders may not be able to 
express their preferences before a specific decision is made. The 
studies considering this type of uncertainty are not in the scope of 
this review paper. 

New markets and new production technologies: Recently, there 
has been an increased attention to utilizing forest biomass for 
bioenergy and biomaterials production. Investment grants, carbon 
and energy taxes, green certificate schemes, conversion technol- 
ogies, and availability and quality of biomass resources may not be 
known with certainty [72]. For example, in designing and planning 
a biomass power plant, it may be hard to estimate the long term 
availability, quality and cost of biomass. Alternatively, market 
demand for biofuels may not be lucid from the beginning. 

In general, uncertainty can be dealt with at the source, or it can 
be dealt with during the process of decision making. When 
uncertainty is ignored, the decision making is based on the 
expected values of stochastic parameters, which may be different 
from their actual values and may lead to non-optimal or infeasible 
results and solutions. Considering uncertainty in decision making 
usually helps companies safeguard against threats on one hand 
and take advantage of the opportunities that higher levels of 
uncertainty would provide on the other hand. It also makes 
decisions robust and mitigates the effect of the variations on the 
optimal solution. 


3.1. Modeling approaches 


In the literature, several approaches were used to incorporate 
uncertainty in the supply chain design and management including 
scenario-based approaches, sensitivity analysis, stochastic optimi- 
zation, and robust optimization. Mula et al. [73] provided an 
overview of the uncertainty in production planning. Based on 
their classification, the supply chain planning and management 
problems with uncertainty are usually solved by conceptual, 
analytical or artificial intelligent based approaches. Some of the 
methods for dealing with uncertainty are out of the scope of this 
paper since they were not employed in forest biomass/bioenergy 
supply chain studies, such as fuzzy sets and fuzzy logic [74], expert 
systems, reinforcement learning, neural networks, genetic algo- 
rithms and multi-agent systems [73]. 

Often times, stochastic models are difficult to solve computa- 
tionally. In many studies, a simple approach is taken by replacing 
the random variables by their expected values. In other words, in 
these studies the developed model is implemented for a single 
future “scenario”. Another simplified method is to generate differ- 
ent scenarios for future and solve the model for each scenario one 
at a time and then combine these different solutions by some 
heuristic rules. Scenario generation itself is a challenging task 
which can be done using historical data, forecasting methods, 
managerial and expert judgment, etc. [75,76]. The scenario-based 
approach provides an extensive what-if analysis which helps in 
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evaluating the outcomes based on different scenarios of the 
stochastic parameters. A set of optimal values attained from this 
method can be used for extracting expected value and standard 
deviation of the objective function and decision variables. How- 
ever, it may not determine the optimal solution for some realiza- 
tion of uncertain parameters or there might be a possibility to have 
a better overall optimal solution. Another issue with this approach 
is that the likelihood of occurrence of each scenario has a great 
effect on the optimal solution, but it cannot be determined with 
certainty. 

The analytical approaches (or stochastic models) are usually 
based on mathematical programming and optimization methods. 
In stochastic modeling, also called recourse programs, it is 
assumed that accurate probabilistic descriptions of the random 
variables such as probability distributions, densities or other 
probability measures are available. Details of stochastic and 
recourse modeling can be found in [75]. In these methods, decision 
variables are divided into two groups called the first stage 
variables (control, here-and-now), which are made before the 
realization of the uncertain parameters, and the second stage 
variables (state, wait-and-see), which are taken after the realiza- 
tion of the uncertain variables. The output of such a model is the 
optimal single first-stage policy and a set of recourse decision rules 
determining which second-stage action should be taken in 
response to each random variable. Often, some decomposition 
methods such as the Bender decomposition are used in the 
recourse modeling to accelerate the achievement of the optimal 
solutions [77]. One of the advantages of developing stochastic 
programming models is in their capability to manage risk asso- 
ciated with the company's performance. 

Another approach for incorporating uncertainty in supply chain 
management and planning problems is called robust optimization. 
This method is used when ranges of random variables are known. 
In this method, a function called “protection function” is formu- 
lated for each constraint containing uncertain parameters. The 
protection function is an optimization problem in itself which 
determines the size of buffer required to immunize a constraint 
against uncertainty based on degree of uncertainty of the para- 
meters and degree of the desired protection. After formulation of 
the protection function and transforming it to an optimization 
model, its linear version of the dual model is obtained and 
embedded into the original model. The resulted model, called 
the robust counterpart, is larger but still linear and can be solved 
by LP solver packages [78]. 

In summary, sensitivity analysis is very useful in determining 
the sensitive uncertain parameters, those parameters that even a 
small change in them would affect the solution. Then, models 
considering variations only in the sensitive parameters, instead of 
all the uncertain parameters, could be developed. The scenario 
analysis technique, which is a highly effective and easy to under- 
stand approach, can be used to deal with uncertainty. In this 
method, a series of cases that represent the possible input 
parameters are considered independently. The solution (output) 
for each scenario is determined separately and compared with the 
solution of other scenarios based on a performance measure to 
make the required decision. When the set of scenarios is large, a 
Monte Carlo simulation can help in generating additional statis- 
tical information for the solution, such as average, frequency, 
bounds, etc. Stochastic programming also considers a set of 
scenarios, however, unlike scenario analysis, it considers all of 
the scenarios at once with each scenario at a certain probability 
and finds a solution that is feasible for all scenarios. When the 
range of uncertain parameter is known, not its distribution, then 
the robust optimization method can be employed. Fixed minimum 
and maximum ranges for uncertain inputs are usually used in this 
method. Robust optimization models can be solved efficiently as 


the size of the problem grows slightly when uncertainty is added. 
Moreover, it is suitable for cases where scenarios or distribution 
functions cannot be defined and only a range for uncertain 
parameter can be acquired. 


3.2. Previous studies 


There are few studies which considered uncertainty in the 
forest biomass supply chain and all are reviewed here. One study 
which used sensitivity analysis in the forest biomass supply chain 
was conducted by Kim et al. [79]. They developed a mixed integer 
linear programming optimization model for the supply chain 
design of bio-gasoline and biodiesel production from six forestry 
resources (logging residuals, thinnings, prunings, inter-cropped 
grasses, and chips/shavings). The first set of conversion plants 
could be from a set of candidate sites with four capacity options to 
convert biomass to three product types: bio-oil, char and fuel gas. 
These intermediate products could be used either as local fuel 
sources or as feedstock to produce final products (gasoline and 
biodiesel) at the second conversion plants, which could be from a 
set of candidate sites with four capacity options. There were 
several possible markets for the final products with certain 
maximum demands. The objective of the model was to maximize 
the overall profit by determining the number, location, and size of 
the conversion plants, biomass supply locations, the logistics and 
the amount of materials to be transported between the various 
nodes of the designed network, while satisfying the demand 
constraints. The considered case study was based on an industrial 
database related to a case in the Southeastern United States. The 
authors evaluated the trade-off between centralized and distrib- 
uted network designs. They also run their model for different 
demands (100% to 90%, 75% and 60%) to evaluate the effect of 
changes in demand on the optimal network design. The results 
showed the total profit for the distributed system was higher than 
that for the centralized design at 100% demand. When the demand 
decreased, the profit difference between the two systems reduced 
as well. In this study, although demand uncertainty was evaluated 
through a very simple sensitivity analysis method, it showed the 
importance of considering uncertainty as it affected the optimal 
result. 

Kim et al. [80] performed a global sensitivity analysis and two- 
stage stochastic programming on a biofuel supply chain model and 
evaluated the effect of uncertainty in different parameters on the 
final result using Monte-Carlo simulation. The uncertain para- 
meters included availability and acquisition cost of each biomass 
type, annualized capital cost of two conversion processing plants, 
cost of transporting intermediate products, biomass and final 
products, value of each intermediate product at the first conver- 
sion processing site, yield of intermediate products from biomass 
at the first conversion processing, yield of final products from 
intermediate products at the second conversion processing, max- 
imum demand, operating costs for two conversion processing 
plants, and the price of each final product. The authors concluded 
that the most important uncertain parameters affecting the profit 
were the price of the final product, the conversion yield ratios of 
the two conversion processes, maximum demand and biomass 
availability. They then generated 33 scenarios from changing these 
five most important uncertain parameters by +20% and developed 
an optimal model to maximize the expected profit from all these 
scenarios plus the expected value scenario. The first stage decision 
variables were the size and location of the processing plants and 
the second stage recourse decision variables were flows of biomass 
and product within plants and markets. They implemented the 
robustness analysis and Monte Carlo global sensitivity analysis to 
compare the performance of the multiple scenario design with the 
single scenario design. It was concluded the impact of variations in 
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stochastic parameters on the optimal solution was mitigated in the 
model that optimized multiple scenarios. In this paper, uncertain 
parameters were only changed within some range rather than 
using probability distributions. Moreover, it only modeled the 
supply chain in one time-step and ignored the possibility of having 
correlation between random parameters. 

Strategic decisions regarding choosing investment options for heat 
savings and decreasing energy imports or increasing energy exports in 
pulp mills under market uncertainty was studied in [81]. The objective 
function of the developed model was to maximize the expected net 
present value of the investments. The decision variables were related 
to investment in heat savings and energy conversion technologies as 
well as distribution of the obtained heat from different energy 
conversion technologies. The uncertainty was considered in future 
energy prices and policy instruments through a scenario tree of five 
different combinations of several emission reduction policies, and 
electricity, lignin and bark prices. The original optimization model was 
a mixed-binary linear programming model and uncertainty was 
included using a multistage stochastic programming model. Each 
stage contained 5 years and the total planning horizon was assumed 
to be 30 years. In [82], Svensson and Berntsson introduced a 
methodology for making decisions about investment in new energy 
technologies in industrial plants. The main focus of the paper was to 
include uncertainty in the energy market. The authors used stochastic 
programming and scenario analysis to take into account the uncer- 
tainty in CO, charge, and crude oil and electricity prices. They then 
applied the methodology to a pulp mill and presented examples of 
possible future investments to show the usefulness of the proposed 
approach. 

Chen and Fan [83] developed a mixed integer stochastic 
programming model to incorporate uncertainty in strategic plan- 
ning of bioenergy supply chain systems. They considered bioetha- 
nol production, feedstock procurement, and fuel delivery in an 
integrated model to minimize costs. The raw material included 
both agricultural and forest biomass. A two-stage stochastic 
programming model was developed and applied to a case study 
in California. Uncertainty was considered in available feedstock 
supply and fuel demand. A Lagrangian relaxation based decom- 
position solution algorithm, called progressive hedging method, 
was used to reduce the computational effort needed to solve the 
stochastic model. This method decomposed a stochastic problem 
across scenarios by partitioning the original problem into manage- 
able sub-problems. The first stage decision variables were refinery 
and terminal locations and sizes and the second stage decision 
variables were feedstock procurement and transportation, ethanol 
production and transportation. In the baseline scenario, four 
discrete demand scenarios with equal probabilities were consid- 
ered and the results showed that the stochastic programming 
provided lower expected and more reliable (within smaller range) 
total system cost than deterministic solutions. In the second 


Table 5 


Summary of studies on optimization of forest bioenergy supply chain with uncertainty. 


Author Year Uncertain parameter 
Kim et al. [79] 2010 Demand 
Kim et al. [80] 2011 Price of the final product, the conversion yield ratios 


of the two conversion processes, maximum demand 
and biomass availability 

Svensson et al. [81] 2011 Emission reduction policies, electricity, lignin 
and bark prices 

Svensson et al. [82] 2011 COz2 charge, fuel and electricity prices 

Chen and Fan [83] 2012 Feedstock supply and demand 


Tay et al. [84] 2012 Raw material supply and product demand 


Method 
Sensitivity analysis 


Two-stage stochastic programming 


scenario, uncertainty in biomass supply was considered using 
ten scenarios of feedstock with equal probabilities and a fixed 
demand. The authors stated that the optimal solution of the 
bioethanol supply chain was not sensitive to the uncertainty in 
supply. This study did not integrate different sources of uncertainty 
within a single framework. It also ignored uncertainty in production 
yields and prices. In addition, the model had single annual time- 
step and did not consider multiple shorter time variations. 

One of the by-products of the Kraft process in the pulp and 
paper industry is black liquor, which can also be used for energy 
production. Tay et al. [84] considered uncertainties in raw material 
supply and product demand in integrated biorefineries using the 
robust optimization method. The model was mixed integer non- 
linear programming with one time step. Different scenarios with 
associated probabilities were considered for supply and product 
demand. The results identified the optimum capacity of each 
process technology and its corresponding amount of biomass, 
intermediate and final products. 

Table 5 summarizes the previous studies which incorporated 
uncertainty in the forest biomass supply chain management and 
design. 

The studies reviewed here indicated the benefits and advantages of 
incorporating uncertainty in the decision making over developing 
deterministic models only.Generally, there are several challenges in 
designing and solving a model with uncertainty which have prevented 
its wide spread applications. First, these mathematical models are 
usually complex. Second, considering uncertainty in the models 
usually requires more time and solving them requires more computa- 
tional effort compared to the deterministic cases. When uncertainty is 
considered in several parameters, the size of the problem might 
become too large to be solved. Another challenge is the difficulty in 
quantifying uncertainty in mathematical terms. In case of stochastic 
models, estimating the probability distributions of each random 
parameter and the probability of occurrence of each scenario are not 
straightforward and may have great impact on the final solution. If 
some random variables are correlated, the estimation of their prob- 
ability is even more difficult. Statistical data analysis and forecasting 
methods which may be based on the historical data along with 
managerial judgment could be used in estimating such probabilities. 
Generating adequate finite scenarios for modeling an uncertain 
environment with infinite possible outcome (in case of continues 
variables) is another challenge. 


4. Conclusions 


Energy generation from forest fuel has several advantages, 
most importantly in terms of mitigating adverse environmental 
impacts of fossil fuels. It helps societies diversify their energy 
sources by providing local energy for communities and 


Case study 


Biofuel industry in 

Southeastern United States 

Biofuel industry in Southeastern United 
States 


Multi-stage stochastic programming Pulp industry in Sweden 


Multi-stage stochastic programming Pulp industry in Sweden 
Two-stage stochastic programming 
with decomposition 
Robust optimization 


Biofuel industry in California 


Kraft pulp and paper industry 
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potentially sell bioenergy products in the energy market. Com- 
mercialization of this energy source, however, is still under 
consideration. Characteristics of the forest biomass supply chain 
for energy production, such as unpredicted raw material quality 
and costly transportation, make it among expensive energy 
sources in the market. Optimization of forest biomass supply 
chains could help in commercializing this sustainable energy 
source by reducing costs. 

Mathematical programming and optimization techniques were 
used in the design and management of forest biomass supply chains. 
These models were used to provide the optimum solution for 
decisions related to the network design, technology choice, plant size 
and location, storage location, mix of products and raw materials, 
logistics options, supply areas, and material flows. The objective 
functions included profit/cost, CO emissions, travel time, etc. Both 
deterministic and stochastic mathematical programming models were 
developed for forest biomass supply chains. 

As environmental and social impacts of products and processes 
are becoming more important than the pure economic value, 
future studies using deterministic modeling should incorporate 
environmental and social objectives in addition to the economic 
objective into the models and perform multi-objective optimiza- 
tion. Considering non-economic objectives into the models would 
help deal with important issues, such as carbon emissions, land 
use, communities' interests, job creation, governmental policies, 
ecological impacts of removing residues from forest areas, and 
recreational aspects of forests, and provide possible trade-offs 
between different objectives. 

Deterministic models are necessary and helpful but not enough for 
capturing all aspects of forest biomass supply chains. Most of the 
parameters and events in forest biomass supply chains, such as wood 
quality, market situation, prices, and yields, are uncertain. Uncertainty 
makes this industry volatile and the long term planning difficult. Even 
in the short term, variations in raw material quality can impact the 
performance and efficiency of forest bioenergy production. Despite the 
interest and much effort in extending deterministic models, only a 
number of previous studies considered uncertainty in optimizing 
forest bioenergy supply chains. Understanding the problem character- 
istics, gathering sufficient data and choosing the appropriate metho- 
dology are important in developing stochastic models in forest 
biomass supply chains. The appropriate methods to be implemented 
depend on the characteristics of the problem, data availability and 
form, and type of uncertainties. 

In the literature, sensitivity analysis and stochastic program- 
ming approach were used to deal with uncertainties in supply, 
demand, prices, conversion yields, carbon tax and emission 
reduction policies in forest biomass supply chains. None of the 
previous studies dealt with uncertainty in biomass quality, such 
as moisture content and heating value, and its effect on the 
produced energy and its costs, especially transportation costs. It 
is worthy to integrate different forest supply chains due to their 
dependency and model the whole value chain considering 
harvest areas, sawmills, pulp mills, and bioenergy conversion 
facilities, and evaluate the effect of uncertainty in one chain on 
the supply, production and demand of the chains. This, however, 
would make the modeling and solution approach more compli- 
cated. Correlation between uncertain parameters is another 
interesting topic which has been overlooked in the design and 
management of forest biomass supply chains. 
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